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I The integration procedure based on the the generalized Darboux transform is suggested for 

the Ishimori magnet model. Exact solutions are constructed for the model of background of 
D ■ spiral structures. The possibility of phase transition in the system is hypothesized. 

in _ 

PACS: 52.35.Sb 

C/} \ It is well known that the phenological approach suggested by Landau and Lifshitz 

in the theory of ferromagnetism is based on the idea that the evolution of weakly ex- 
cited states of a spin system in the long- wavelength limit can be described in terms of a 
constant-length magnetization vector (magnetic moment density) and is characterized by 
a certain effective field [1]. With allowance for the fact that the greatest contribution to 
the process comes from the exchange interaction between crystals atoms, this made it pos- 
sible, in particular, to obtain nonlinear equation for the one-dimensional case (isotropic 
Heisenberg ferromagnet model), which was subsequently solved by the inverse scattering 
transform method (ISM) [2, 3]: exact solitonic solutions were obtained, excitation spec- 
trum was described, and infinite set of conservation laws was found. Further progress in 
this direction was achieved in work [4], where the Lax representation was applied to the 
anisotropic Landau-Lifshitz model (and, thereby, it was proved to belong to the class of 
completely integrable models), and in [5], where its exact solutions were found by the 
"dressing" method. Noteworthy is also work [6], which was devoted to the construction 
of the integrable deformations of Heisenberg model. 

The situation is much more complicated in the two-dimensional case. The corre- 
sponding nonstationary equation for a Landau-Lifshitz ferromagnet, both isotropic and 
anisotropic, proves to be nonintegrable (see, e.g. [7]) and has exact solutions only for 
rather specific cases. At the same time, there is much evidence, both experimental [8] and 
obtained by numerical simulation [9], of the existence of stable localized two-dimensional 
excitations with a finite energy. For a two-dimensional system, the spectrum of these exci- 
tations becomes more diversified; in particular, new nontrivial topological objects appear 
in the spectrum. 

The model suggested by Ishimori in [10] is presently the major tool for the phenomeno- 
logical description of ferromagnets of dimensionality (2+1). This is, primarily, due to the 
key property of this model: it is completely integrable and allows the use of ISM and 
<9-dressing procedure to construct a rather broad class of solutions (vortices (lumps), ra- 
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tional exponential solutions, instantons, etc.) on the trivial background [11-14], while, 
by using a nonstandard Darboux transform [15], a physically interesting solution can be 
obtained in the form of a vortex circulating on circumference with a constant angular 
velocity. 

However, it is worthy of note that the exact two-dimensional solutions were found at 
the expense of necessity of introducing nonlocal (along with the exchange) spin interac- 
tions into the model. The physical mechanism of this interaction is as yet unclear. In 
this connection, it should be emphasized that the standard Heisenberg (exchange) inter- 
action mechanism in the so-called Schwinger-boson mean-field theory has become more 
understandable only a relatively short time ago (see, e.g., [8] and references cited therein). 
This gives grounds to hope that the nature of nonlocality that provides a broad spectrum 
of reasonable solutions and the corresponding observed physical objects will be clarified 
in the future and that the Ishimori model is quite realistic (within the framework of the 
adopted macroscopic approach). 

In this work, a new and rather effective integration procedure is proposed for the 
Ishimori ferromagnet model. It opens the way for constructing exact solutions, including 
the ones on the nontrivial background. As to the ISM (or the <9-dressing method), this 
may encounter considerable technical difficulties. 

The Ishimori magnet model is given by 

S t = S A (S xx + ^Syy) + UyS X + U X Sy, (la) 

u xx - a 2 u yy = -2a 2 S(S x AS y ), (lb) 

where S(x, y,t) = (S±, S 2 , S 3 ) is the three-dimensional magnetization vector, |S| = 1, u = 
u(x,y,t) is the auxiliary scalar real field, and the parameter a 2 equals ±1. The case 
a 2 = 1 will be called Ishimori-I magnet model (MI-I), and a 2 = — 1 will be called MI-II. 

Note that in the static limit (S = S(x,y)) and u = const, the MI-I model transforms 
into the model of two-dimensional isotropic Heisenberg ferromagnet (elliptic version of 
the nonlinear 0(3) - sigma-model), which was integrated in [16-17] by using the ISM with 
boundary conditions of the spiral-structure type. 

A characteristic feature of model (1) is the presence of topological charge 

-i poc poo 

Qt = — \ / S(S X A S y )dxdy, (2) 

J — oo J — oo 

which is conserved in the course of system evolution (integral of motion) and represents 
the mapping of a unit 2-sphere onto the 2-sphere: S 2 — > S 2 . This mapping is known 
to be characterized be the homotopic group ^(S* 2 ) = Z, where Z is the integer group, 
signifying that Qt must be integer. According to Eq.(lb), the scalar function u = u(x,t) 
is related to the topological charge production density 2 . 

The set of Eqs.(l) is integrated using the following associated linear system: 

tf„ = -S* x , (3a) 
a 

2 Although, strictly speaking, the function u has no direct physical meaning, the functions u y and 

u x , related to each other by Eq.(lb), can likely be interpreted as "frictions coefficients" inducing forced 

precession of the magnetization vector S along the the x and y axes, respectively. 
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% = -2iSV xx + Q* x , (3b) 

where Q = u y I + a 3 u x S + iaS y S — iS x , \I/ = fy(x,y,t) G Mat(2,C), S = Yli=i Si a ii °% are 
the standard Pauli matrices, and / is a unit 2x2 matrix. From its definition, the matrix 
S has the following properties: S = S*, S 2 — I, det S — — 1 and Sp S = (the symbol 
(*) denotes Hermitian conjugation). 

In what follows, we restrict ourselves to the MI-II case (a = i); the situation with 
MI-I can be analyzed in a similar manner. 

We will solve Eq.(l) by the method of generalized Darboux matrix transform 3 . We 
also demand that system (3) be covariant (U — > U, ^ — > \&) about the transformation of 
the form (U = S) 4 

* = n(*,* 1 )*?, (4) 

where *i = ^i(x,y,t) is a certain nongenerate bare solution to the system of Eqs.(3), 
and Q(^f, ^i) = Q(x,y,t) G Mat(2,C) is a functional defined on the pair set of matrix 
functions 5 . One then obtains from Eq.(3a) two dressing relations: 

u = nv 1 u^ 1 n- 1 , u = in y (n x )- 1 . (5) 

Hence, setting z = x + iy, z = x — iy, d z = (1/2) (d x — id y ), d s = (l/2)(d x + id y ), 
and = ^^U^i, we obtain the first equation for the matrix il (det (I ± U) = 

det(I±WW)) = 0): 

(J-H/(°))(fi- 1 ) 2 = 0. (6) 
Note that, due to the symmetry relation U = —a 2 Ucr 2 , the following involutions are valid: 

qr = a 2 ^fa 2 , Q = a 2 Qcx 2 , 0, = cr 2 Cla 2 , (7) 
which means, in particular, that Eq.(6) can be rewritten in the "conjugate" form 

(J + ^KQ^^O. (8) 

To obtain the second equation for the matrix fl, one should use Eq.(3b). This equation, 
however, is rather cumbersome. Taking into account that Eqs.(5) and (6) can easily be 
expressed in terms of function the following consideration can best be carried out for 
this function. 

Using the identity Q + UQ = -2i(I + U)(ul + U) s , one obtain from Eq.(3b): 

(/ + U){% + 2W XX + 2i(u- z I + U^ x } = 0. (9) 

Let us transform this equation. Making allowance for the relation q> yz - = —U^ X — iU^ xS 
that follows from (3b), one has 

3 In [18], the approach based on the Darboux transform was applied to a system that is gauge equivalent 
to (3), and the dressing relations were obtained thereafter. 

4 In the literature, the scalar variant of this transformation is sometimes called Moutard transform 
[19]. 

5 It follows from Eq.(3a) that the matrix solutions \& and are related to each other by the nonlinear 
relationship ^("J^)" 1 = #ij,(#i x ) _1 . 
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(/ + U){% + 2i$ xx + 2iu s * x - 2^ yz - 2iUV x ,} = 0. (10) 

After multiplying both sides of this equation by U on the left and adding together the 
resulting expression and Eq.(10), one finds: 

(/ + U){% + 2i^> zz + 2W- Z - Z + Aiu z ^ z } = 0. (11) 

Thus, after applying transformation (4) together with the covariance requirement, one 
obtains two (and two analogous conjugate) equations for the function 

(/ - U)* s = 0, (12) 
(/ + U){% + 2M ZZ + 2M gs + Aiu z ^ z } = 0. (13) 
Next, taking into account the identities S(S X A Sj,) = (l/(2i)) Sp (UU x U y ) and 

Sp (UU x U y ) = Sp (UU x U y ) + 2iA(lndet + 2iSp {[U, C4]* -1 * x + 

(14) 

+[U, Uy]*- 1 ^} + 4Sp {Ufi-^y, 

where A is the two-dimensional Laplace operator, and requiring that Eq.(lb) be covariant, 
we rewrite it as 

A(u - u - 2 In det *) = 2 Sp{[U, U x ]^>- 1 ^ x + [U, Uy^'^y - 2iU[^- 1 ^ y , AT 1 **]}. (15) 
The relevant expression for the dressed topological charge has the form 

Qt = Qt H / / dx dy A In det 

J — CO J -co 

+2-/ / dxdj/Sp {[^,^]*- 1 ^+[^,^*- 1 i y -2iC/[^- 1 ^,*- 1 * a; ]}. (16) 

«/ —CO -V —CO 

Introducing the notation [1] = C/[l] = U,..., u[l] = u, . . . ,Q T [1] = Qt, ■ ■ 

one can readily obtain from Eqs.(5), (15), and (16) upon N-fold dressing of the starting 
bare solution U = [/« = U[0] 6 : 

(N-l \ /N-l \ _1 



f + OO p+oo 

u[N] = u + 2 ^lndet + / / dx'dy' G(x - x',y - y')x 

j = l J —oo J —oo 



(18) 



6 These formulas can also be written in terms of the matrix functionals ^ j), where ^i, are 

certain bare solutions to the set of Eqs.(3) with S = S^K 
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N 



X 



J'=l 



g T [iV] = g T + 



47T 



^2SpAj(x', y', t), 

dxdy^A In det * [j] + 



j'=i 



(19) 



1 

47T 



+00 



N 



-1 : 



where A,-(x, y, f) = [f/[j - l],U x [j - * x [j] + [U[j - 1], C/,[j - l]j v^]- 1 *^] - 

2it7|j - 1] *D'] _1 * V |>'], ^[j]" 1 ^^] , = (1/2tt) ln(a; 2 + y 2 ) is the Green's function 

of the Laplace equation. 

Turning back to the simple dressing, we note the Eqs.(12) and (13) defines the whole 
collection of the solutions to model (1) (in reflectionless ( in ISM terms) section of the 
problem). According to Eqs.(12) and (13), four cases is possible: 1). ^f z = and the 
braced expression in Eq.(13) is also zero; 2). = and the expression in braces is 
nonzero; 3). the situation opposite to (2); and 4). ty z 7^ and the expression in Eq.(13) 
is also nonzero (i.e., the columns of this matrix belong to the kernels of the degenerate 
I — U I + U transformations, respectively). 

In this letter, we restrict ourselves only to the first case. Then, 



tt* = 0, * t + 2i* zz = 0. 
This system has the well-known polynomial solutions (\& 
*n, *i2 = -*2i)([10], [13]): 



(20) 



= 1,2, * 
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Ni 

E E 

j=0 m+2n=j 



a 3 



1 



Mi 



mini 



-it)", * 21 (M) 



E E 

j=0 m+2n=j 



mini 



(21) 

where N\, M\ are natural numbers and M\ — N± — 1, a,j and bj are the complex numbers; 
and the first summation is over all possible combinations of the numbers m, n > such 
that m + 2n = j. 

The bare solution to the set of Eqs.(l) is taken in the form of vector function = 
(0, sinc&W, cos^ 1 )), where = 5 t + a x + Poy + Jo, a o, Po, 7o? ^0 £ K. are parameters), 
i.e., in the form of a two-dimensional spiral structure with Qt = = (according to 



Eq.(2)). To determine the function u(x, y, t) = w '(x, y, t) one should substitute this 
vector into Eqs.(la) and (lb). The requirement for the compatibility of the resulting two 
linear equations gives, after integration, 



u 



(i) 



g£\y + ?°z) + fg?W) + ^x(s'),t) ds>, 
a J a 



(22) 



where g^ an d 9? are arbitrary functions. The function is constant on the charac- 
teristic y + (f3o/a>o)x = const and s is its parameter. Therefore, the explicit expression 
for u is determined by Eq.(18) with N — 1 and by Eqs.(21) and (22). 
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Dressing relation (17) gives (S + — Si + iS 2 ) 

COS$ (1) (|^n| 2 - |^2l| 2 ) - i Sin (§21^11 - *21*11; 



Ss(x, y, t) 



|*ll| 2 + |^2l| 2 



(23) 



5 ( 2cos$ (1) ^2i^n + isin$ (1) (^ 2 1 + & 21 ) 
b + {x, y, t) — -= — -= — . 

*n 2 + * r 12 



'21 1 

Setting 5o = and N± — 1, we obtain the simplest static anti (vortex) (one-lump) solution 
on the background of (also) static spiral structure: 



S 3 (x,y,t) 



cos$«[l 


a - |aizj 


2 


6 | 


2 ] + i sin $W [6 


(a - |aiz) - b (a - \a x z)\ 




<2g — ^CliZ] 


2 + 


6 | 


2 



(24) 



2b cos (a - \a x z) + i sin [6q + (a - \a\zf\ 



S + (x,y,t) - , , , 

|a - |ai-z| 2 + |o | 



The calculations in Eq.(16) show that — * oo; the divergence arises after the integration 
of the first two terms in the braces. For 5 = a = (3q = 70 = 0, solution (24) transforms 
into a static (anti) vortex (on the trivial background) with the topological charge Qt = — 1 
(see also [13]). 

For 6 ^ and N t = 2, i.e. for vj/ n = a - (a 1 /2)z + (a 2 /2)[(l/4)^ 2 - it] and 
^21 = b — (l/2)6i2;, formulas (23) describe the dynamic (anti)two-vortex (two-lump) state 
on the background of (also) dynamic spiral structure with Q T — > 00, and it transforms to 
the state with Qt = —2 if the parameters entering turn to zero. 

Clearly, the intermediate types of solutions are also quite realistic and of interest. 
Among these are a static vortex on the background of the dynamic spiral structure and 
a dynamic vortex on the background of the static spiral structure. 

One the basis of these results, the hypothesis can be put forward that the struc- 
tural second-order spiral-vortex — > vortex phase transition (analogous to the Kosterlitz- 
Thouless transition) is possible in the system of interest. This can occur if the parameters 
Sq, cto, A) and 70 are functions of time, i.e., functionals of an external nonstationary and 
spatially uniform magnetic field 7 . Then the fact that the parameter turns to zero means 
that there is a certain critical field ("Curie point" or, more precisely, Lifshitz point) that 
corresponds to the phase transition point. This hypothesis is confirmed by the experimen- 
tal fact that the spiral (modulated, incommensurate) structure in a magnetic field can 
convert into the commensurate structure corresponding to a paramagnet with magnetic 
moments mainly oriented along the external field 8 [20]. We also emphasize that this phase 

7 One can show that the addition of the term S A H(t) to the right-hand side of Eq.(la), where H(t) is 
the external magnetic field, merely renormalizes the magnetization vector; i.e. this term can be eliminated 
by the appropriate transformation through making the bare solution dependent on magnetic field. 

8 Clearly, the theoretical justification of this hypothesis should rest on the consideration of the order 
parameter of system and on the analysis of the Ginzburg-Landau functional in the vicinity of the critical 
point. However, although the Hamiltonian of system (1) is known (generally speaking, it is obtained in 
[21] for a modified MI model), the relevant calculations become overly cumbersome even for the simplest 
solutions and, thus, are beyond the scope of this article. 
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transition should be accompanied by a change in a symmetry and topological properties 
of the system. 

Another series of solutions to model (1) can be found if the solution to the system of 
Eqs.(20) is sought in the form 

/oo 
B 11;21 (p)e- 2ip2t+pz dp. (25) 
-oo 

Here Bu^i are the functional parameters. In particular, by setting -811,21 (p) = cn^iS(p — 
Pn,2i), where S(.) is the Dirac delta function and cn^i G C and pu^i G K. are parameters, 
one gets using Eq.(17) (cn ;2 i 7^ and the symbol c.c. stands for the complex conjugation): 

S 3 (x,y,t) = 



COS$ (1) [ Cn 


2 e 2|pn|x _ 


C21 


2 e 


2\P21 


\ x ] — i sin 


<J> (1) [ciiC2ie 2i ^?i +p 2i) t+2 ^ 11+P21 ) 2 - c.c] 




Cll 


2g2pna; _|_ 


C21 


2g2p2iz 



(26) 

S+fo y, *) = 

2cnC2i cos^e 2 ^"-^!^ 11 ^ 21 * + isin$( 1 )(c? 1 e 4ip?1 * +2pi1 * + c 2 .^- 4 ^* 2 *' 21 *) 

| Cii |2 e 2 Pll a: + | C2l |2 e 2p 21 z 

Therefore, a exponential and nonsingular solution is found on the background of spiral 
structure. At — > 0, the solution is finite if pn, p 2 i > 0; in this case the component 
S3 evolves only along the x variable. 

More complicated solutions of this type can be found if the functionals -811,21 are taken 
as linear combinations of delta functions. 

Of interest is to compare the results obtained in this work with the results of work [22] , 
where the MI-II model was proved to be gauge equivalent to the known hydrodynamic 
Davey-Stuartson-II system (that describes the evolution of nearly monochromatic small- 
amplitude wavepacket at the surface of a small depth fluid). This means that the Lax pair 
can be transformed from one system to another by a certain gauge transformation, which, 
in turn, allows the formulas relating the solutions for these systems to be derived. It is 
significant that the initial boundary-value problems with specified classes should posses 
similar equivalence; in [22], the class of rapidly decreasing Cauchy data was assumed in 
both cases. Clearly, there is no gauge equivalence in the considered case of spiral (and, 
hence, nondecreasing) structures, while the solution constructed on their background 
cannot be derived from the solution to the Davey-Stuartson-II model. 

Note in conclusion that the approach developed in this work can easily be extended 
to a series of Myrzakulov magnet models [23,24], which are modifications of the Ishimori 
model; for them, the first Lax-pair equation either is close or coincides with Eq. (3a) and 
the main modifications concern the functional Q in Eq. (3b). 

I am grateful to A.B.Borisov for attention to the work and S.A.Zykov and A.V.Shirokov 
for assistance. 
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